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Abstract. The paper is devoted to the study of the motion of one-dimensional rigid bodies 
during a free fall in a quasi-Newtonian hyperviscous fluid at low Reynolds number. We show 
the existence of a steady solution and furnish sufficient conditions on the geometry of the body 
in order to get purely translational motions. Such conditions are based on a generalized version 
of the so-called Reciprocal Theorem for fluids. 



1. Introduction 

The study of the free fall of slender bodies in liquids is an essential issue in many problems 
of practical interest, such as the design of composite materials or the analytical technique of 
separation of macromolecules by electrophoresis (see [3] for a very interesting and rich review 
on the subject). Typical experiments show that homogeneous bodies satisfying some symmetry 
conditions, when dropped in a quiescent viscous liquid, will eventually reach a steady state that 
is purely translational, having the symmetry axis forming an angle with respect to the gravity g, 
called tilt angle, that depends on the material geometry of the body and on the physical properties 
of the liquid. 

If the geometry of the body is such that one of the dimensions dramatically prevails on the 
other two, the assumption that the body is one-dimensional is a reasonable simplification which 
can give satisfactory results. However, a one-dimensional body is "too thin" to interact with a 
classical Newtonian incompressible fluid in 3D (it has null capacity, see [9]). Hence, we propose to 
study the problem of the free fall of a slender body in a regularized model for Newtonian fluids, 
introduced by Fried and Gurtin [2] in 2006, where higher-order derivatives are considered in the 
constitutive prescription of the Cauchy stress tensor. 

The Navier-Stokes equation for incompressible fluids reads 

Pttt + Pi u ■ V)« + Vp - pAu = pb, 
at 

where p is the pressure field, u is the divergence-free velocity field, p > is the constant and 
homogeneous mass density, p > is the dynamic viscosity and pb is a volumetric force density. 
The hyperviscous regularization consists of adding a term proportional to AAm to the equation. 
For this modified equation, 

du 

p— + piu ■ V)u + Vp - pAu + (AAu = pb, 

at 

where £ > is the hyperviscosity, the existence and uniqueness of regular solutions have been 
established. In a series of papers [2, 8, 6, 4] a purely mechanical explanation of the hyperviscous 
term has been proposed, and different contributions to £ associated with dissipation functionals 
are introduced and analyzed. Here we assign to £ a geometric role, by introducing the effective 
thickness L > of the lower-dimensional objects, and setting £ = pL 2 , so that the hyperviscous 
flow equation becomes 

du 



p— + p{u ■ V)u + Vp - pA (u - L 2 Au) = pb. 
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In the experiments on the free fall of rigid bodies in viscous fluids the Reynolds number is often 
very small, so that the inertia of the liquid can be neglected and one can linearize the flow equation 
[9, 1]. However, even after that approximation the problem does not become fully linear, since 
there remains a nonlinear coupling between the flow and the rigid body motion. 

In the present paper we will study the steady fall of one-dimensional rigid bodies in a hyper- 
viscous fluid at low Reynolds number. In Section 2 we give the mathematical formulation of the 
problem, and in Section 3 we study the forces acting on the body in the case of a hyperviscous 
fluid. In Section 4 we show the existence of steady solutions, and in Section 5 we prove the Recip- 
rocal Theorem in the case of a linearized hyperviscous liquid surrounding a one-dimensional body, 
and study some properties of the so-called resistance tensors. Finally, Section 6 contains sufficient 
conditions on the geometry of a homogeneous body in order to get purely translational solutions. 



2. Formulation of the free fall problem 

The free fall problem is characterized by the fact that the rigid body is immersed and dropped 
from rest in an otherwise quiescent fluid and gravity is the only external force acting on the system. 
We represent a rigid body as a connected, bounded, closed subset S of M 3 which is a finite union of 
images of [0, 1] through a C 1 -diffcomorphism. It is convenient to write the problem in a co-moving 
frame, with origin at the center of mass c(t) of S. Denoting by y the position of a point in the 
original inertial frame, and by x its position in the co-moving frame, we know that, at any time 
t > 0, 

x = QT(t)(y-c(t)), 

where Q(t) is an orthogonal linear transformation for any t > 0, with Q(0) = 1. If the velocity of 
the center of mass and the spin of the rigid body in the inertial frame are denoted by T](t) and 
17 (i), respectively, so that 

v(t) = V (t) + fl{t) x(y- c(t)) 

is the velocity, in that frame, of any point y belonging to the rigid body, then their expression in 
the co-moving frame is given by 

£(f) := Q J (t)r](t) and := Q T (0^)> 

respectively, and the rigid velocity field v is transformed into 

U(t) := £(t)+u(t) x x, 

where x denotes the coordinates of y in the co-moving frame. 

The gravitational acceleration vector g (constant in the inertial frame) is represented, in the 
co-moving frame, by G(t) := Q T (t)g, which is easily seen to satisfy the ordinary differential 
equation 

— = Gxu;. (1) 

As customary when studying flows past rigid bodies, the velocity field u that we consider is the 
so called disturbance field, which is the difference between the actual flow and the flow at infinity, 
both seen in the co- moving frame. Since the flow at infinity is — U (that is minus the extension to 
all of the fluid of the motion of the immersed object), the representation of the fluid flow in the 
co-moving frame is given by u — U. 

The continuity and flow equations for an incompressible (disturbance) velocity field u{x,t) and 
pressure field p(x, t) defined on (K 3 \ S) x [0, +oo), become then 

divM = 0, (2) 

( du \ 
p ( — + [(it -U) ■ V]it + uj x uj = div T(u, p) + pG, (3) 

where T(u,p) denotes the Cauchy stress tensor. Notice that, thanks to its frame indifference 
properties, T retains the same functional dependence on the velocity field seen in both the inertial 
frame and in the co-moving one. 
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The disturbance field u satisfies also the decay condition 

lim u(x,t) = 0, (4) 

\x\— too 

and the adherence to the rigid body, given by 

u(x, t) = U(x, t) onSx[0,+co). (5) 

To properly account for Archimedean forces we introduce the effective mass of the body as 
given by 

m e = to — m c , 

that is, the difference between the real mass of the object and the complementary mass m c of 
a portion of fluid occupying the real volume of the object. Even if for a one-dimensional body 
the complementary mass m c should vanish, in view of the interpretation of such bodies as rep- 
resentations for real three-dimensional objects we allow for any value < m c < m, suggested by 
the physical properties of the interaction between the body and the fluid. Then the equations of 
motion for the rigid body in the co-moving frame read 

d^ 

to— + muj x £ = m e G + / (it,p), (6) 
J^+(jx(Ju) = -m c rxG + t(u,p), (7) 

where J is the inertia tensor of S, r is the position of the centroid' 1 -* in the co-moving frame, and 
/ and t are the total hydrodynamic force and torque, respectively, exerted on the rigid body as a 
consequence of the fluid flow u and pressure p. The proper definitions of / and t, as well as the 
physical properties of the fluid encoded in the Cauchy stress tensor T, are discussed in Section 3. 

The whole set of equations (l)-(7) represents the differential problem associated with the free 
fall of a rigid object S in an incompressible fluid. It is convenient to consider it in the non- 
dimensional form: by choosing suitable reference speed W and length d, we can switch to non- 
dimensional quantities according to 

X t[l U pd TO 
X — > — , t — > , U — > , p — > , TO — > , 

d ' pd 2 ' W y fj,W pd 3 ' 

divw = 0, (8) 

— + Re{[(u - U) ■ V}u + uj x u} = divT(it,p) + G, (9) 
lim u(x,t) = 0, (10) 

|a3|^-oo 

u(x,t) = U{x,t) on S x [0,+oo), (II) 

to^ + Re(mu) x £) = m e G + f(u,p), (12) 
dt 

du) 

J— + Re\u3 x (Jw) = -m c r xG + t(u,p), (13) 
dt 

^ = Re(G x W ), (14) 
with initial conditions 

u(x,0) = £(0)=u>(0) = 0, G(0)=g, (15) 

where Re = pWd/p is the Reynolds number and every quantity has to be understood as non- 
dimensional. Moreover, a suitable choice of the reference quantities ensures that |G(£)| = 1 for 
any t > 0. 



obtaining 



W Notice that the centroid of a rigid body coincides with its center of mass when the body has a uniform mass 
density; in the latter case r = 0. 
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The low-Reynolds-number approximation of the differential problem, which is also a lineariza- 
tion of the equation for the flow, is obtained by neglecting the terms proportional to Re in equa- 
tions (9), (12), and (13), whereas equation (14) remains unchanged, since it represents a geometric 
constraint which holds for any non- vanishing value of Re. Moreover, the steady version of the prob- 
lem is achieved by assuming that all the quantities do not depend on time, hence neglecting all the 
time derivatives. In that case, keeping into account also the initial conditions (15), equation (9) 
becomes 

divT(it,p) +g = 0, (16) 

and equation (11) writes 

u(x) = U(x) = £ + lj x x on E. 
Finally, equations (12)— (14) become 

m e g + f{u,p) = 0, -m c r x g + t(u,p) = 0, gxu = 0, 

respectively. 

3. The viscous force acting on a slender body 
Given r > 0, we introduce the (closed) r-neighborhood of the slender body E by setting 



V r (Y) := I x G R 3 : inf \x - c| < r 

Then we define the total hydrodynamic force, due to the fluid velocity and pressure field (u,p), 
acting on E as 

f(u,p) := lim / T(u,p)n, (17) 



dV r (T,) 



where n denotes the unit outer normal to dV r (Y,). Notice that, thanks to the regularity of E, the 
r-neighborhood V r (T,) has a Lipschitz boundary for any r sufficiently small. 

Proposition 3.1. The limit in (17) is well-defined. 

Proof. We consider a ball Br centered at the origin and with radius R, which contains V r (Yi) 
for some r > 0. According to equation (16), the term divT balances the gravity, so that it is 
represented by a measure whose singular part is concentrated on E (we can see this by noting 
that the mass density per unit volume must diverge on E to give a non-zero weight to a body 
with vanishing volume). Denote by A ac and A s the absolutely continuous and singular parts of 
the measure div T, respectively. It follows that the support of A s is contained in E and that 
divT = A ac in R 3 \ E. Then we have, by applying Lebesgue's theorem, 

lim / divT= lim / X ac = / A ac = / divT, (18) 

JB R \V r (£) JB R \V r (Z) Jb r \£ JB r \S 

since A ac G L 1 (R 3 ;IR 3 ). Then, by the Divergence theorem, 

/=lim / Tn= / In - [ divT, (19) 

where n is always the outer normal. Since the right-hand side of (19) is independent of r, the 
left-hand side is well defined. □ 

In a similar fashion, we define the total hydrodynamic torque acting on E, due to the fluid 
velocity and pressure field (u,p), as 



t{u,p) :~ lim / xx T(u,p)n. 



ldV r (T,) 

It is important to stress the fact that, if A s were absent, i.e. if divT were an L 1 -function, then 
the integral over Br \ E of div T in (19) would be equal to its integral over all of Br and / and t 
would simply vanish. 
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The constitutive theory for non-simple fluids leading to a hyperviscous flow equation has been 
developed in [2, 8, 6]. It offers a number of possible choices for the terms to be included in T, in 
addition to those of Newtonian fluids. Here we make the following, somewhat minimal, choice: 

T{u,p) := -pl+ (Vm + Vm t -£ 2 VAu) . (20) 

In this way we obtain a fluid which is quasi-Newtonian, while being able to adhere to lower- 
dimensional objects. The only new parameter £ is given by L/d, hence it is non-dimensional 
and strictly positive; as explained in [5], it can be assigned the geometric meaning of an effective 
thickness of the slender body S. It is straightforward to check that T, as defined in (20), enjoys 
the standard symmetry and frame indifference properties and satisfies a dissipation inequality. In 
particular, one has 

div T = -Vp + Am - £ 2 AAu. 

We summarize the problem of the steady free fall of a one-dimensional body £ at low Reynolds 
number, as the following: find (u,p) G £f 2 (R 3 ;R 3 ) x if" 1 ^ 3 ; M) and £,u,g G R 3 with \g\ = 1, 
such that 

divu = 0, (21) 
Vp - Am + ^ 2 AAtt = g on R 3 \ S, (22) 
u(x) = £ + w x x on S, (23) 
m e g = -/, (24) 
m c r x g = t, (25) 
g xu = 0. (26) 

The set H k is the usual Sobolev space of functions with square-intcgrable weak derivatives up 
to the order k, and the set H~ k is the topological dual of H k ; the strong decay condition (4) is 
encoded as an integrability condition in the assumption u G if 2 (R 3 ;R 3 ). Notice that, although 
equation (22) is linear, the full problem is nonlinear. 



4. Steady free fall at low Reynolds number 

In this section we prove the existence of a solution for the the differential problem (21)-(26). 
We begin by introducing some auxiliary problems, which are well-posed by virtue of the following 
result. 

Lemma 4.1. Given £, ui G R 3 , there exists a unique solution (h,p) G -ff 2 (R 3 ;R 3 ) x iJ _1 (R 3 ;R) 
of the problem 

divh, = in R 3 , 

Vp - Ah + £ 2 AAh = in M 3 , (27) 
h = £ + ui x x on S. 

Moreover, such a solution is of class C°° on I 3 \ E. 

Proof. Since 7? 2 (R 3 ;R 3 ) embeds in a space of Holder-continuous functions, the subset 

H := {v G ff 2 (R 3 ;R 3 ) : div v = and v = £ + w x x on £} 

is well-defined, closed and convex. The velocity field h can be found by minimizing on H the 
functional 

^H-^l [ (2|SymVu| 2 + £ 2 \Av\ 2 ) , 

where SymVu := (Vt> + Vi> T )/2. Being T a strictly convex functional, h is unique. Then, the 
pressure field p can be recovered as the Lagrange multiplier of the divergence-free constraint. 

The same technique used to study the regularity of the solution for the classical Stokes equation 
can be applied here to prove that the solution is of class C°° on R 3 \ S. □ 
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Consider now the solutions (h {i) ,p^) and (iJ (i) ,P«) (i = 1,2,3) of 

'divh (i) =0 inR 3 , 

VpW - A/i (i) + £ 2 AAh {i} = in R 3 , (28) 

% h (i) =ei on S, 
and 

'div.Ff (l) =0 inR 3 , 

VPW - Ai? (l) + £ 2 AAH {i) = in R 3 , (29) 

= e,- x x on E. 



We will show that the combinations 

3 

u = J2lZih (i) +UiHW), p = Y^iP^+uJiP^]+g- X , (30) 

i=l i=l 

for a suitable choice of the vectors £ and w, solve the steady free fall problem. First we need to 
introduce four matrices, which will play an even more important role in Sections 5 and 6. 

Definition 4.2. The matrices K, S, C, and B arc defined in Cartesian components by 

Kji := - lim / T(h {l} ,p {t) )n ■ e h (31) 



Sh := - lim / T(H U \ P ij) )n ■ e t , 

r ^JdV r (T.) 



(32) 



C^:=-Um/ x x T(h u \p (j) )n ■ e u (33) 

r ^°JdV r (Z) 

Bji := - lim / xx J{H {l \P {l) )n ■ e h (34) 

r ^°JdV r (S) 

where n is the outer normal to V r (T,). Following [1], they are called resistance tensors, and in 
particular K is the translation tensor, B the rotation tensor, and S and C the coupling tensors. 
Moreover, we denote by A the 6x6 matrix 

'K S N 
C B y 

We postpone to Theorems 5.3 and 5.4, in the next section, the proof of a fundamental property: 
the matrices K, B and A are symmetric and positive definite. 
Although an energetic argument of Brenner [7, Section 5-2] is usually adopted in this case, we 
will prefer to give a direct proof. 

Now we can prove the main theorem of the section. 

Theorem 4.3 (Existence Theorem). The differential problem (21)-(26) admits a solution (u,p, £, oj, 

Proof. It is straightforward to check that the fields u and p defined by (30) satisfy equations (21), 
(22), and (23). Equation (26) implies that u> = Xg for some A S R, and equations (24) and 
(25) reduce to the following algebraic system in the six scalar unknowns £, A, and g (recall that 
l9l = l): 

K£ + ASg = m e g 
C£ + ABg = —m c r x g. 

It is now clear that the steady free fall problem admits a solution if and only if (35) admits a 
solution, and the latter fact is related to the properties of the matrix 

K S 
C B 

Since K is non singular, the first equation of (35) becomes 

C=K- 1 (m e9 -AS ff ) 
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and one can eliminate £ in the second equation of (35). Since A is non-singular, the linear 
transformation 

Fg := (CK _1 S - B)- 1 (m e CK _1 fl + m c r x g) 
is well-defined and non-singular, and we can write (35) as 

/ t = K-\m e g + >Sg) 

\ F 9 = Xg. (36) 

Being F a 3 x 3 real matrix, it has at least one real eigenvalue. Such an eigenvalue A, the associated 
unit eigenvector g and £ calculated as in the first equation of (36), together with the fields u and 
p introduced in (30), furnish a solution for equations (21)-(26). □ 

5. An analysis of the resistance tensors 

In (31)-(34) we introduced the four resistance tensors K, S, C and B. In view of the conditions 
on E assumed in the auxiliary problems (28)-(29), we can give an equivalent characterization. 

Proposition 5.1. The resistance tensors are such that 

Ku = - lim / T(/i w , p« )n-h {3 \ 

Sji = -lim / J{H {3 \P^)n-h {t \ 

Cji = - lim / T(/i w , )n ■ H (l) , 

Bji = -lim / T(H {l \P^)n- H u \ 

r ->°JdV r (£) 

Proof. Since h^' is smooth and = on S, one has, as r — > 0, 

- eiW^ev^ := sup \h®(x) - e;| -> 0. 

Hence, considering for instance the translation tensor K, it follows that 
f T{h {l \p^)n-h u) - f T(h«,pW)n-« 



< \\h& - e^mm I T(fcW,pW)n 

JdVUS) 



>dV r (E) 

as r — > 0. The proof of the other three formulae is similar. □ 

In the remainder of the section we will prove that K and B are symmetric and that C T = S. 
We first need a fundamental property of steady incompressible flows at low Reynolds number, 
the so-called Reciprocal Theorem (see [7, Sec. 3-5]), which, roughly speaking, states a reciprocity 
property between two solutions of the same equation, independently of the boundary conditions. 
The validity of the theorem, which is quite trivial for ordinary fluids, is not so obvious in the 
present case of hyperviscous fluids, since the lack of further boundary conditions and the higher 
order of the differential operator can break such a reciprocity. However, the theorem can be 
recovered for the particular case of one-dimensional bodies. 

Theorem 5.2 (Reciprocal Theorem). Given two solutions (u\,pi) and (1x2,^2) in the space 
H 2 (R 3 ;R 3 ) x ir-^R 3 ;^) of the equations 

divtt = inM 3 , divT(w,p) = in M 3 \ E, (37) 

where T{u,p) is defined as in (20), we have 

lim / T(ux,pi)n ■ u 2 = lim / T(u 2 ,P2)n ■ u\. 
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Proof. Consider a large ball Br containing V r (T,) and apply Gauss-Green formula to the domain 
B R \Vr(Z): then 

/ T(u 1 ,p 1 )n-u 2 = - T(ui,pi) ■ Vm 2 - / T(ui,pi)n • u 2 , 

JdV r (T,) JB R \V r (Z) JdB n 

where the normal in the left-hand side is exterior to V r (T,) and we kept into account that 
divT(u2,P2) = on M 3 \ E. The last surface integral on OBr vanishes as R —> +00, since 
any solution to the hyperviscous Stokes' problem decays as l/\x\ (see Appendix A), hence T 
decays as 1 / 1 | 2 . 

Now take the first term of the right-hand side and use the constitutive prescription (20): 



- / T(«i,pi) • Vm 2 

JB R \V r (T,) 

= -/ (Vwi + Vit{) • Vm 2 +t 2 / VA«i • Vw 2 

= - I (Vu 2 + VmJ) • Vtti +£ 2 VA«i • Vm 2 - 

JB R \V r {S) JB R \V r {T.) 

Consider the last term; since the gradient and Laplace operators commute, by Green's second 
identity it follows that 



VA«i • Vit 2 = / Vtti • VAm 2 

B R \V r (£) JB R \VrCS) 



- / [Vit 2 • (VVui)n - Vmi • (Wu 2 )n] 

- / [Vit 2 ■ (Wui)n - Vui • (Wit 2 )n] . 

J OB R 



Now we claim that the surface integrals vanish as r — > and R — > +00. Indeed, take for instance 
the term 

/ Vit 2 • (Wtti)n 
and apply again Gauss-Green formula inside V r (Y,). Then 

/ Vit 2 • (VViti)n = / [VVtti • VVtt 2 + V« 2 • AV«i]. 

JdV r {T,) Jv T (T,) 

Since u u u 2 G i? 2 (K 3 ;K 3 ), then VVwi • VV« 2 6 L^R 3 ) and 

lim / VV«i • VV119 = 0. 

In a similar way, since Vit 2 G ff^R 3 ; Mat 3x3 ) and AViti £ H^ 1 (R 3 ; Mat 3x3 ), then V« 2 AVmi e 
L 1 ^ 3 ) and 



lim / 



Vtt 2 • AV«i = 0. 

'V r (E) 

Regarding the surface integral on 9-Br, it vanishes as R — > +00 since it 1,1*2 decay as l/|a?| (see 
Appendix A). □ 

By combining Proposition 5.1 with Theorem 5.2 one can immediately prove the main result of 
the section: 

Theorem 5.3. The resistance tensors are such that 

K T = K, B T = B, C T = S. 
In particular, the matrix A is symmetric. 

Now we are in a position to give a proof of the positive definiteness of the tensors K, B and A. 
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Theorem 5.4. The 6x6 matrix A is positive definite. As a consequence, also K and B are positive 
definite. 

Proof. For £ = (£1,62,6) and 10 = (w ls w 2 , w 3 ) set 

i=l i=l 

where (ft ,p^) and (H^\P^) are the solutions of the auxiliary problems (28)-(29). Using 
Proposition 5.1, the Reciprocal Theorem and the linearity of T one can check that 



CO 



£ ■ K£ + 2co ■ C£ + w • Bo; = - lim / T(u,p)n • u 



Now argue as in the proof of the Reciprocal Theorem. Take a large ball Br containing V r (S) and 
apply Gauss-Green formula to the domain Br \ V r (Yj) to obtain 

— / J{u,p)n-u = I T(u,p) • V«+ / J{u,p)n-u, 

JdV r {£) JB R \V r (E) JdB R 

where the normal in the left-hand side is exterior to V r (T,) and we kept into account that 
divT(u,p) = on M 3 \ E. The last surface integral on 8Br vanishes as R — > +00 since the 
solution decays as l/|a;|, hence T decays as l/|a;| 2 . 

Now consider the first term of the right-hand side and use the constitutive prescription (20): 

/ T(m,p)-Vm= / (Vm + Vm t ) • Vu - I 2 / VAm-Vu 

JB B \V r (E) JBh\K(S) •/B H \V r (E) 

|Vm| 2 + / Vu T -Vu-l 2 / VA« • Vu. (38) 

B R \V r (S) JB R \V r (T,) JB R \V r (T,) 

We deal with the second integral, taking into account that divit = 0: 

/ Vu T • Vu = - / u ■ (Vu T n) + / u- (Vu T n). 

JB R \V r (X) JdVr(E) JdB R 

The last integral on 8Br vanishes as R — > 00 for the usual asymptotic behavior at infinity. By 
applying the Gauss-Green theorem on the other integral it follows that 

- / u ■ (Vtt T n) = - / Vtt • Vtt T - / u ■ div(Vu T ) 

JdV r (T,) JV r (S) iv r (S) 

and both terms vanish as r — > since the first integrand is in L 1 and in the last integral one has 
div u = 0. 

Now we deal with the last integral of (38). Since the gradient and Laplace operators commute, 
we have 

-i 2 [ va« -vu = i 2 [ 1 vv«| 2 

JB R \V r (£) JB R \Vr(S) 



[ Vu ■ (VVu)n -£ 2 [ Vu- {VVu)i 



idVriT.) JOB 

Following the last part of the proof of the Reciprocal Theorem, it can be proved that the surface 
integrals vanish as r — > and R — > +00. Summarizing, 



CO 



A(*]=/ \Vu\ 2 +e 2 \VVu\ 2 , 

JB R \V r (T,) JB R \V r (T.) 



hence A is (strictly) positive definite. □ 
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6. Translational solutions for bodies with symmetries 

The free fall of a one-dimensional body in a hyperviscous fluid at low Reynolds number is char- 
acterized by 21 independent coefficients: 12 coefficients for the tensors K and B and 9 coefficients 
for the coupling tensor C. However, material symmetries of the body can significantly reduce 
such a number. Moreover, the symmetries induce some restrictions on the form of the resistance 
tensors. We are specifically interested in symmetries which induces purely translational motions 
of the body (that is, with oj = 0). We will now study some particular symmetries. 

If the body is invariant under a change of frame given by an orthogonal matrix Q, then also the 
solutions of the auxiliary problems (28) do not change; on the contrary, the solutions 
of (29) undergo a sign change if det Q = — 1, due to the presence of the vector product in the 
boundary condition. Hence one can prove that 

K = Q T KQ, B = Q T BQ, C = (det Q)Q T CQ. (39) 

6.1. Bodies with a plane of symmetry. We say that the body E has x 2 a; 3 as a plane of 
material symmetry, if the density function p of the body satisfies 

p(-Xx, X2,x 3 ) = p(xi,x 2 ,x 3 ) for every {x\, x 2 , £3) € r3 . 

In particular, a homogeneous body has a plane of material symmetry if, and only if, it is symmetric 
with respect to that plane. 

Proposition 6.1. Assume that E has x 2 x 3 as a plane of material symmetry. Then the resistance 
tensors have the form 

K = 

Proof. Since E is invariant under the orthogonal transformation given by 

/-I 0\ 
Q = 1 , dctQ = -1. 
\0 ij 

then formulae (39) yield 

K12 = K13 = B i2 = B13 = 0, 
Cn = C 22 = C33 = C 2 3 = C32 =0. □ 

Now consider the system (36) which solves the problem of the steady free fall, in the case 
when the body E has a plane of material symmetry, say £ 2 x 3 . Suppose moreover that E is 
homogeneous, so that the center of mass and the centroid coincide, hence r = 0. In such a case 
the second equation of (36) becomes 

TO e (CK- 1 C T - B)- 1 CK" 1 g = Xg. (40) 

Since uj = Xg, we get a translational solution whenever A = 0. Being (CK _1 C T — B) and K positive 
definite matrices, (40) has a solution A = if, and only if, det C = 0. In the case of a body with 
a plane of material symmetry, indeed, the latter condition is satisfied and it is easy to check that 
an eigenvector of C, say Uo, corresponding to the eigenvalue A = lies in the plane X2X3. Hence 
one has 

m e (CK- 1 C T - B)- 1 CK" 1 g = <=> g = Ku ; 

by the form of K given in Proposition 6.1, also the vector Kuq lies in the plane X2X3. We can 
summarize the latter result in the following theorem: 

Theorem 6.2. Assume that E has X2X3 as a plane of material symmetry. Then there exist 
an orientation of the body, lying in the same plane of symmetry, which gives rise to a purely 
translational solution. 



Now it is quite easy to study the class of bodies with two orthogonal planes of symmetry: 
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il 



Corollary 6.3. If the body has two orthogonal planes of symmetry, say x\x% and X2X3, then the 
free fall along the x^-direction gives rise to a purely translational motion. 

Proof. Since the body is invariant under the orthogonal matrices 
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-1 
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it is easy to check that K and B are diagonal, and C has the form 





^0 


C12 


H 


C21 







^ 






Hence Uq = (0,0, a) is an eigenvector of C corresponding to the null eigenvalue, and the motion 
with orientation given by 

furnishes a purely translational solution. □ 

6.2. Helicoidally symmetric bodies. Now we study bodies which are invariant under the action 
of a rotation of angle 9 G [0; 27r[ around the 21-axis, which is represented by the orthogonal matrix 




Following [7], we say that a (one-dimensional) body S is helicoidally symmetric if there exists a 
co-moving frame such that 

Rg£ = £ for some 9 ^ 0, 7r, 

that is, if it is invariant under a discrete group of co-axial rotations of order strictly greater than 2. 
For instance, a homogeneous body composed by three concurrent edges of a regular tetrahedron 
is helicoidally symmetric with 9 = 2n/3. 

Proposition 6.4. Assume that E is helicoidally symmetric around x\. Then K and B are diagonal 
with K22 = K33 and B22 = B33, and C is of the form 













H 





C22 


C23 






— C23 


C33 



Proof. Let us employ formulae (39) with Q = Rg, keeping into account that det Rg = 1. For the 
matrix C we get the conditions 

C12 = C12 cos 9 + C13 sin 6>, 
C13 = C13 cos 9 — C12 sin (9, 
C21 = C 2 i cos 9 + C31 sin 6>, 

C 22 = C 22 cos 2 9 + (C 23 + C 32 ) cos 6> sin + C 33 sin 2 6>, 
C23 = C23 cos 2 9 + (C33 - C 22 ) cos 6> sin - C 32 sin 2 6>, 
C31 = C31 cos 9 — C21 sin0, 

C32 = C 32 cos 2 9 + (C33 - C 22 ) cos 9 sin 9 - C 23 sin 2 6», 
C33 = C 33 cos 2 9 - (C23 + C 32 ) cos6»sin0 + C 22 sin 2 6», 

which in turn imply that C12 = C21 = C13 = C31 = and C23 + C32 = 0, since 9 ^ 0, ir. Being K 
and B symmetric, we have the further conditions K23 = B23 = and K22 = K33, B22 = B33. □ 



12 



G. G. GIUSTERI, A. MARZOCCHI, AND A. MUSESTI 



6.3. Helicoidally symmetric bodies with fore-aft symmetry. A remarkable situation is the 
case of a homogeneous one-dimensional helicoidally symmetric body with fore-aft symmetry, that 
is, a body which is both helicoidally symmetric around an axis, and has a plane of symmetry 
orthogonal to that axis. 

Without loss of generality, let us assume that a one-dimensional body X is helicoidally symmetric 
around x\ and has X2X3 as a plane of symmetry. Since the coupling tensor C has to satisfy both 
Proposition 6.1 and Proposition 6.4 at the same time, it follows that 

C = 0. 

Assuming that £ be homogeneous, so that r = 0, the system (36) merely becomes 

U = m e K- 1 g, 
[0 = Xg, 

hence A = for any direction g. Then u> = and for any given orientation the body falls with a 
purely translational velocity given by £ = m e K~ 1 g. 

Appendix A. Green's function for Stokes flow 

The basic tool used to construct solutions to the Stokes problem is the so-called Stokeslct, that 
is the Green's function for the Stokes operator in R 3 . In this Appendix^ 2 - 1 we want to compute 
the expression of the Stokeslct in the case of our hyperviscous fluid, identified by the operator 

A := £ 2 AA - A. 

We first need a Green's function g solution of the fourth-order elliptic equation 

£ 2 AAg - Ag = S(x - x'). 
Using the Fourier transform, we easily obtain 

^ r e ik-(x-x') 
9( X - X ') = / |.i 2 /, 2 |. l2 , ^ dk. 



(27T)3 J r3 |fc|2(^|fc|2 + 1) 

We choose a basis for the momentum space in such a way that x — x' is along the fc 3 -direction, 
set R = I a; — x'\, switch to polar coordinates (k, 8, </>), and use the calculus of residues to obtain 

- = (0^2 1 \ x WTW d{cos9)dk 

2 f +co sinfci? 

dk 



(2ir) 2 PRJ k(k 2 + l/P) 
1 



Im 



Hence, the Green's function is 

g(x - x') 



(2ir) 2 £ 2 R 
1 

{2tt) 2 £ 2 R 



1 

4ir\x — x' 



+00 ikR 

dk 



k{k 2 + l/l 2 ) 



2 n£ 2 e~' 



1 — exp 



f. 



(41) 



Notice that, in the limit £ — > 0, (41) reduces to the fundamental solution 

gi (x - x') = — -;— (42) 

for the Laplace operator. Moreover, g is well defined for any x G R 3 , at variance with the classical 
expression g±, which is singular at the origin. 

We now proceed to construct the hyperviscous Stokeslet, that is a pressure field and a velocity 
field £ satisfying 

div C = 0, (43) 



results of the Appendix are based on [5]. 
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Vp c +AC = hS(x), (44) 

with h G K 3 and A = £ 2 AA — A. Let <j> satisfy A<j> = 8(x — x 1 ); then, since A commutes with V, 
a solution for (43)-(44) is given by 

p c = -Ad, £ = h4> + W. 

The scalar field d entering this solution is chosen to satisfy the constraint (43) and turns out to 
have the explicit form 

= (-A)" 1 ^ • V<t>) = gi * (h ■ V<f>), 

where g\ is as defined in (42) and * denotes the usual convolution product. Now, exploiting 
the properties of the convolution and the operator A, and using the Green's function g given by 
equation (41), we find that 

-(.9i * A{h ■ Vff)) = -(.9i * div(A(gh))) = - div( 5 i * A(gh)) = h ■ V 9l 

and, denoting by / the Fourier transform of the function /, 

91* (K- Vfl) = J i(h ■ k)g l9 e^ x dk 

i(h-k)e lkx 



(2tt) 3 
1 

(2^ 
h ■ x 

4ir 2 £ 2 \x\ 
—h ■ x 
4tt 2 £ 2 \x\ 
hx rl 



|fe| 4 (^ 2 |fc| 2 + l 



■dk 



ik\x | cos 6 



x k(k 2 + i/e 2 ) 

ho ° sin(fc|a|r) 
o k{k 2 + l/P) 



-d(cos 6)dk 



dkdr 



8n\x\ 
h ■ x 

~8n\x\ 



Hence the Stokcslct is given by 



p c {x) 



t\ (l-e- 




2t i-i 


2£ 2 


e i 

N 


+ w 


h ■ x 




4ir\x\ 3 


1 



(It 



- 1 



C(x) = 



h 



8irx 



, „ _M 21 _M 2i 2 
1 — 2e r — - — -e J - - — - 



e~— - 1 



(ft, • x)x 

87r|a;| 3 



l + 2e™ 



— p * 



6^ 2 



- 1 



We also define the hyperviscous Oseen tensor Z as, using Cartesian components, 



Zij(x) :— 



87TSC 



, „ i-i 21 m 
1 - 2er~ - t—^~ 
\x\ 



2t 2 

w 



* - 1 

l + 2e"V 4 



U _w 6^ 2 / 

n e T + m e 



- 1 



8tt|cc| 3 

whereby it follows that C(x) = Z(cc)ft. The Stokeslet allows us to obtain an integral representation 
for the solution of 

Vp - A(u - £ 2 Au) = pb, 
with vanishing condition at infinity, in the form of a convolution: 



u(x) ■= p I Z(x — x')b{x')dx' . 
Jm 



(45) 



In particular, whenever b has compact support, such as in the case of the gravity force acting on 
a bounded body E, the solution u behaves as l/\x\ for \x\ — > oo. 
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